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Abstract 

In this paper we apply the equation of state (EOS) of QCD at finite chemical potential and zero 
temperature proposed in H. S. Zong and W. M. Sun [Int. J. Mod. Phys. A23, 3591 (2008)] to 
the study of properties of quark star. This EOS contains only one adjustable parameter m/j which 
sets the scale of chiral symmetry breaking (in our calculation we have chosen two values of m/5: 
m/3 = 244 MeV and m^ = 239 MeV, which is fitted from the value of /jr and determined by 
e^e~ annihilation experiment, respectively). From this EOS a model of quark star is established by 
applying the Tolman-Oppenheimer-Volkoff equation under two conditions: with the P{fi = 0) term 
and without the P{fi = 0) term. Our results show clearly that the P(^ = 0) term is an important 
quantity in the study of quark star. A comparison between our model and other models of quark 
star is made. In particular, we have compared our results with the most recent observational 
data measured using Shapiro delay reported in P.B. Demorest et al. [Nature (London) 467, 1081 
(2010)]. 
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I. INTRODUCTION 



Since Gell-Mann and Zweig [l, 2| proposed the conception of quark, many authors have 
used this idea in the study of astrophysical problems js-S]. Just as the discovery of neutron 
led to the idea of neutron star, the conception of quark naturally stimulated many physicists 
and astrophysicists to suggest that quark star maybe exists as well js-llj. In the later 
decades the research about the nature of quark star was developed. As is well known, the 
equation of state (EOS) of matter plays a key role in the study of the structure of compact 
stars. However, at present it is very difficult to obtain the EOS of compact matter from the 
first principle of Quantum Chromodynamics (QCD) in all ranges of density and temperature. 
Therefore one has to adopt various QCD inspired models to study the properties of dense 
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ects. For example: the MIT bag model 
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-|l9l|. and the perturbative QCD model |20l-|23| . Although all of these models can refiect 



15| . the Nambu-Jona-Lasinio (NJL)model 



some basic features of QCD, none of them can provide a complete and consistent framework 
to work in all ranges of density and temperature. The MIT bag model can provide the 
quark confinement mechanism, the NJL model can describe chiral symmetry breaking of 
QCD, while the perturbative QCD model works well at high energy scale due to asymptotic 
freedom. Present astrophysical observations and laboratory experimental data cannot enable 
one to make a clear choice among these models. So it is worthwhile to develop other models 
with better QCD foundation to study compact matter. Recently, in Ref. [24] the authors 
derived a model-independent formula for calculating the EOS of QCD at finite chemical 
potential and zero temperature with the aid of functional integral formalism and obtained 
an explicit analytic expression for the EOS by applying the rainbow-ladder approximation 
of the Dyson- Schwinger approach and the model quark propagator proposed in Ref. [25j]. In 
this paper we shall apply this EOS to the study of structure of quark stars by considering two 
situations: with the P{fi = 0) term (the pressure at zero chemical potential and temperature) 
and without the P(/i = 0) term. 

This paper is organized as follows. In Sec. HIl we first give a brief introduction to the 
model-independent formula for the EOS of QCD at finite chemical potential and zero tem- 
perature derived in Ref. [24] and give the explicit analytic expression of the EOS obtained 
using this formula with the rainbow-ladder approximation of the Dyson- Schwinger approach 



and the model quark propagator proposed in Ref. 



25| . From these results the properties of 



quark star, such as mass-radius and mass-energy- density relations are obtained by solving 
the Tolman-Oppenheimer-Volkoff (TOV) equation under two conditions: with the P(/i = 0) 
term and without the P(/i = 0) term. A comparison between our model and other models 
of quark star is made. In Sec. IIIIl we offer a discussion of our study. 



II. STRUCTURE OF QUARK MATTER 



(i) EOS of quark matter 



It is well-known that in equilibrium statistical field theory 26| , once the partition function 
(or equivalently, EOS) of a system is obtained, then all the thermodynamical variables can 
be determined. The calculation of the EOS of QCD at finite chemical potential is a contem- 
porary focus. However, due to dynamical chiral symmetry breaking and confinement, it is 
very difficult to get some valuable and reliable results about EOS of QCD from perturbation 
theory at finite chemical potential. So people have developed various nonperturbative meth- 
ods for studying strong interaction physics at finite density and zero temperature |27H32|. 



In Ref. 



24| ] the authors derived a model-independent formula for calculating the EOS of 



QCD at finite chemical potential and zero temperature with the aid of functional integral 
formalism. The expression of the EOS derived there reads 



p(/x) = p(/i = 0) + r djj'piij') 

Jo 
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where N^. and Nf denote the number of quark colors and quark flavors, respectively; G[^{p) 
and p(yu) = —NcNf J -^^tr{G[fj,']{p)'y4} are the dressed quark propagator and the quark- 
number density at finite /i, respectively. From Eq. (1) it can be seen that in order to 
obtain the EOS at finite chemical potential, a key issue is to know the behaviors of the 
dressed quark propagator at finite /i. At present it is impossible to derive the exact quark 
propagator at finite fi from first principles of QCD. In Ref. 33j a new approximation method 
for obtaining the fermion propagator at finite fi was proposed by using the rainbow-ladder 
approximation of the Dyson- Schwinger approach and then the validity of this method was 
discussed in detail in Refs. 34j-l37l]. According to Ref. 33|, in order to obtain the dressed 
quark propagator at finite /i, one should specify the dressed quark propagator at /i = in 
advance. In this work we choose the following model quark propagator proposed in Ref. 



25|: 



with = 1, 5(p2 



(2) 



where m^) is a parameter connected to a dynamically gener- 



ated quark mass. Then, according to the result in Ref. 
at finite /x can be written as 

1 



^7 ■ pA^p^) + B(p'^) 



the dressed quark propagator 



(3) 



where p = {p^p^ + i^). Using the above G'[/i](p), the authors in Ref. |2J] obtained an explicit 
analytic expression for the EOS by means of the formula ([1]). The quark- number density 
obtained there is 



P(/^) 
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and the pressure -P(/i) is 
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and P(/i = 0) (the pressure at zero chemical potential) is calculated from the Cornwall- 



Jackiw-Tomboulis (CJT) effective action 38l. |39| 



P{^^ = 0) = N^Nf J ^{trln [Go(p)G-^(p)] - itr[l - Go ^(p)G(p)] 
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where Go{p) is the bare quark propagator. Substituting the model quark propagator 
into Eq. ([8]), one obtains 
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Here, we want to stress that mo is connected to the dynamically generated quark mass. 
When one discusses low energy QCD at finite fi, it is evident that the dynamically generated 
quark mass should be /i dependent. So, in principle, uid should depend on /i. However, 
in this work, according to the approximation adopted in Ref. [2^, m/? is taken to be n 
independent (In fact, the value of m^j can be determined experimentally, for example, from 
the value of vr decay constant /tt or determined by e~^e~ annihilation experiment). 

Before we turn to the study of quark star by means of EOS ([5]), it is necessary to have a 
look at EOS dSD itself. 

First let us see the large /i behavior of -P(/i). For fj. > A^^^mj:,, the two step functions 
in EOS equal 1. Using the explicit expressions of /(/i) and g{fi), one can expand the 
pressure in EOS around yU = 00 as follows: 



P(/x) = N,Nf 
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This equation clearly shows that, as /i — )■ 00, -P(/i) tends to the pressure of a massless free 

f /i^, and the deviation from the free quark gas result can be expanded 



127r' 



quark gas Po(/^) 

in powers of the ratio itld/ ^. Here it is interesting to compare our EOS ([5]) with the cold, 

22| . The pressure density to second order in in the 



perturbative EOS proposed in Ref. 



MS scheme obtained in Ref. 



22| is quoted as follows: 
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where G = Gq - 0.536iV/ + A^/lnA^/,G'o = 10-374 ± 0.13 and A is the renormahzation 
subtraction point. The scale dependence of the strong coupling constant as (A) is taken as 



as(A) 
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where u = ln(AVA|g), /3o = 11 - 2Nf/3, /3i = 51 - 19A^//3 and /^s = 2857 - 5033A^f/9 + 
325A^|/27. For Nf = 3, A,^ = 365 MeV. The only freedom in the model of Ref. [22j is 



the choice of the ratio A//i, which is taken to be 2 in that reference. The perturbative EOS 
(fTTj) is applicable only in the chirally symmetric phase, when the chemical potential fi is 
larger than fi^, the chiral phase transition point. From Eq. ( fTTl) it can be also seen that as 
/i — )■ 00, the perturbative EOS also tends to the free gas result, and the correction terms 
are from perturbative contributions. By comparing Eqs. ([5]) and ( 1111) . it can be seen that 



for moderate values of fi there is a quantitative difference between our EOS ([5]) and the 
perturbative EOS (11 II) : for our EOS the correction terms are positive and the pressure is 
larger than the free quark gas result, whereas for the perturbative EOS the correction terms 
are negative and the pressure is lower than the free quark gas result. Here it also should be 
noted that in our model 171^ sets the scale for dynamical chiral symmetry breaking, therefore 
the correction terms in our EOS arise from nonperturbative effects. 

Second, let us see the small fi behavior of our EOS ([5]). From Eq. (2) it can be seen that 
the obtained quark-number density distribution differs significantly from the Fermi distri- 
bution of a free quark gas. Physically this is a consequence of dynamical chiral symmetry 
breaking in the low energy region. We note that when fi is smaller than a critical value 
/io = |4^/^m£), the quark- number density vanishes identically, and the pressure P{fi) equals 
P{fi = 0) (i.e., the pressure of the vacuum). Namely, /i = /iq is a singularity which sep- 
arates two regions with different quark-number densities. This result agrees qualitatively 



with the general conclusion of Ref. [40]. In that reference, based on a universal argument, 
it is pointed out that the existence of some singularity at the point fi = fio and T = is 
a robust and model-independent prediction. It should be noted that any EOS of QCD at 
finite chemical potential and zero temperature with good QCD foundation should at least 
satisfy this requirement and our EOS does satisfy it. 

From the above analysis it can be seen that in both the large n limit and small fi limit 
our EOS has the correct behaviors required by QCD, which is quite different from other 
EOSs, such as the perturbative EOS (1111) . In addition, our EOS (3) has another merit that 
it has a simple analytic form and contains only one adjustable parameter mo- So it can be 
easily applied in practical calculation of quark stars. Therefore, we expect that our EOS (3) 
can provide a possible new approach for the study of quark stars. 



Just as was shown in Ref. j4l|, the P(/i = 0) term in the EOS is an important quantity, 
since it enters the energy density, which is relevant for integrating the TOV equations. In 
order to see explicitly the influence of the P(/x = 0) term to the properties of quark star, 
in the calculations that follow, we shall distinguish two situations: with the P(/i = 0) term 
and without the P{fi = 0) term. A discussion on the physical meaning of this term will be 
given in the final part of this paper. 

Now let us turn to using EOS (3) to study the properties of quark star. In order to see 
the behavior of the pressure with the variation of /x, we plot the P — /x curve in Fig. [TJ 
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FIG. 1. The pressure as a function of chemical potential, relative to the free quark gas pressure 
Pj^^^ = NcNfH^/{12TT^). In the SQMl, B = (WAMMeV)^ 

Here, as in Ref. [25], we take rrij^ = 244 MeV in our calculation. From Fig. [T] it can be 
seen that the pressure decreases with increasing under condition one (with the P(yU = 0) 
term), whereas it increases with /i increasing under condition two (without the P(/i = 0) 
term). The results of Fraga et. al. [22] and SQMl 42] are also shown in this figure. In Fig. 
Pwe make a comparison of the results in the present paper with those in our previous paper 



4l|, Fraga et. al. [22| and the SQMl models |42|. It can be seen that the SQM EOS is the 



stiffest and the renombagSR model EOS 4l| is the softest. This can be easily understood 
since the interaction between particles will soften the EOS of the system, 
(ii) Property of quark star 

Using the above results, one can obtain the structure of quark star by integrating the 
TOV equation: 

dP{r) _ G {e + P){M + Anr^P) 

(13j 



G {e + P){M + 4:7rr^P) 



dr 



dM{r) 
dr 



r(r - 2GM) 



Anr e. 



(14) 



The calculated mass-radius relations and mass-energy density relations are shown in Figs. 
3-6, respectively. Here, in order to see the sensitivity of the mass-radius relation on the 
parameter in our EOS (3), in the following calculation we will choose two values of 
ni£,: ni£) = 244 MeV and rriD = 239 MeV, which is fitted from the value of j25| and 



determined by e~^e annihilation experiment 



43|, respectively. In order to see clearly the 
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FIG. 2. Quark pressure-energy relation. 
TABLE I. Model parameters and calculational results 



Without P{fi = 0) term 


With P{fi = 0) term 


md{MeV) 


M© 


R(km) 


p{g/cm^) 


md{MeV) 


Me 


R(km) 


p{g/cm^) 


244 


2.3 


24 


6.76 X 10^^ 


244 


0.80 


8.0 


10^6 


239 


1.8 


17 


5.22 X 10^^ 


239 


0.87 


8.4 


9.33 X 10^5 



calculated results in Figs. 3-6, we list them in Table I. 

From Table I it can be seen that, for the case without the P(/i = 0) term, for = 
244 MeV, the maximal mass ~ 2.3 Mq, the radius ~ 24 km and the mass density ~ 
6.76 X 10^^ g/cm^, which is similar to the results calculated by means of the EOS proposed 



in Ref. [30|| without the P{fi = 0) term. The curve of the mass-radius relation is the same 
as that of a neutron star. For m^i = 239 MeV, the maximal mass ~ 1.8 Mq, the radius 
~ 17 km and the mass density ~ 5.22 x 10^^ g/cm^. Here we note that our results for 
m£, = 239 MeV are very close to the most recent observational data reported in Ref. 44 1 
(in that reference the mass of the pulsar J1614-2230 is calculated to be (1.97 ± O.O4)M0 
and its radius is between 11 and 15 km). As was pointed out in that reference, this result 
rules out almost all currently proposed hyperon or boson condensate EOS, whereas quark 
matter can support a star of such a mass only if the quarks are strongly interacting and 
are therefore not "free" quarks. Our results support this observational data. For the case 
with the P{ix = 0) term, for m^ = 244 MeV, the maximal mass ~ 0.8 Mq, the radius 
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FIG. 3. Mass function of radius. 



~ 8 km and the mass density ~ 1.06 x 10^^ g/civ?. For mjj = 239 MeV, the maximal mass 
~ 0.87 Mq, the radius ~ 8.4 km and the mass density ~ 9.33 x lO"*^^ g/cm^. For these two 
values of m^, the curve of the mass-radius relation is the same as that of a typical quark 
star. From the above results it can be seen that for the case without the P(/i = 0) term, the 
structure of the star is sensitive to the value of m^, whereas for the case with the P(/i = 0) 
term, the structure of the star is insensitive to the value oi m^. 

From Figs. 5 and 6 it can be seen that the mass density obtained without the P(/i = 0) 
term is less than the corresponding one obtained with the P(/i = 0) term. The mass density 
obtained with the P(/i = 0) term is larger than the corresponding ones in all other models 
(p ~ lO^^^j^). This is because the EOS for this case is soft and leads to a more compact star. 
Here it is interesting to compare the study in the present paper with that in a previous work 



by some of the same authors 



4l|. In Ref. j4l|, the P(/i = 0) term is identified with -P, 



where B is the vacuum bag constant, and B is taken as a phenomenological input, whereas 
in the present paper the P(/i = 0) term is determined self-consistently in the rainbow-ladder 
approximation of the Dyson-Schwinger approach. 



III. DISCUSSION 



As we have mentioned above, in the present paper we consider two conditions: with the 
P(/x = 0) term and without the P(yU = 0) term (the P(/i = 0) term is calculated using the 

9 



0.4 
0.2 
0.0 
-0.2 



chiral mode2(m^ 


=244MeV) 


3R 




chiral mode2(m^ 


=239MeV) 



FIG. 4. Mass-radius relation. The 3R curve is calculated from the EOS proposed in Ref. [3C| 
without the = 0) term. 









chiral model 






(m^=244MeV) 






SQM1 






SQM3 






Fraga(2p) 






chiral model 




I.I.I 


(m_|=239MeV) 





14.0 14.5 15.0 15.5 16.0 16.5 

Log[p(g/cm')] 

FIG. 5. Mass-energy density relation. 



CJT effective action). This term is not the same as the bag constant B in the MIT bag 
model, because it is positive while B is negative! Furthermore, in the calculation of structure 
of stars this term leads to negative energy density at small chemical potential, which is hard 
to understand from a physical point of view. This is one of our motivations for adopting 
the EOS proposed in Ref. [2J] to study the properties of quark star. In addition, in spite of 
the progress made in the study of quark stars, there are still some authors who suspect the 
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?1G. 6. Mass-energy density relation. The 3R curve is calculated from the EOS proposed in Ref. 
30] without the P(/i = 0) term. 



existence of quark star. In Refs. |45l -l47j. the authors think that there is no quark matter 
core in the neutron star. Ozel |48| argued from their observational data that quark star 
cannot exist, because in term of their observational results the mass > 2.1 ± O.28M0 and 
ihe radius R > 13.8 ± l.Si^m. Of course, other authors do not consent to her arguments 



491], and the discovery of meta-millisecond pulsar is in favor of quark star 42j. Here we also 



note that in a recent paper 4J], the authors argue that their observational data support 
the existence of quark star with strong interactions. Therefore, the present astronomical 
observations do not rule out the existence of quark stars. 



To summarize, in this paper we adopt the EOS of QCD at finite chemical potential and 
zero temperature proposed in Ref. {24] to study the properties of quark star. A comparison 
between our result and those in the literatures is made. Applying this EOS, one obtains the 
structure of quark star. The results obtained under the two conditions (with the P(/i = 0) 
term and without the P{fi = 0) term) are very different and the reason for this is analyzed. 
This shows that the = 0) term plays an important role in the study of properties of 
quark star. 
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